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Abstract 



We show that, for both scalar and spinor QED, the two-loop Euler-Heisenberg 
effective Lagrangian for a constant Euclidean self-dual background has an extremely 
simple closed-form expression in terms of the digamma function. Moreover, the 
scalar and spinor QED effective Lagrangians are very similar to one another. These 
results are dramatic simplifications compared to the results for other backgrounds. 
We apply them to a calculation of the low energy limits of the two-loop massive 
N-photon 'all +' helicity amplitudes. The simplicity of our results can be related to 
the connection between self-duality, helicity and supersymmetry. 



1 Introduction: QED and QCD in constant fields 



The Euler-Heisenberg Lagrangian, one of the earliest results in quantum electrodynamics Q [2|, 
describes the complete one-loop amplitude involving a spinor loop interacting non-perturbatively 
with a constant background electromagnetic field. Euler and Heisenberg found the following 
well-known integral representation for this effective Lagrangian: 
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Here T denotes the (Euclidean) propertime of the loop fermion, and a, b are related to the two 
invariants of the Maxwell field by a 2 — b 2 = B 2 — E 2 ,ab = E • B. The subtractions of the terms 
of zeroeth and second order in a, b corresponds to on-shell renormalization. 

The Euler-Heisenberg Lagrangian contains the information on the low - enery limit of the 
one-loop A^-photon amplitudes for any N. Moreover, it does so in a form which is extremely 
convenient for the study at low energies of nonlinear QED effects [|| such as photon-photon 
scattering photon dispersion [Q, 0], and photon splitting |], 0, fi], Qj. For scalar QED, an 
analogous result was obtained by Schwinger Q: 
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The Lagrangians (|Ll]) and ( |1.2| ) are real for a purely magnetic field, while in the presence of 
an electric field there is an absorptive part, indicating the process of electron-positron (resp. 
scalar-antiscalar) pair creation by the field ||. 

The first radiative corrections to these Lagrangians, describing the effect of an additional photon 



exchange in the loop, were obtained in the seventies by Ritus |1C, 11, 12]. Using the exact 
propagators in a constant field found by Fock fl3| and Schwinger H, and a proper-time cutoff 

(2) 

as the UV regulator, Ritus obtained the corresponding two-loop effective Lagrangians C scal spin 
in terms of certain two-parameter integrals. Similar two-parameter integral representations for 

(2) 



seal spin wnere obtained later by other authors, using either proper-time []14|, |15[ or dimensional 



regularisation |16|, 17 1. Unfortunately, all of these double parameter integral representations are 



quite complicated, so that it is much more difficult to study the weak- and strong-field expansions 
at two-loops than at one-loop. Even more difficult becomes the analysis of the imaginary part of 



the effective action [18]. This is true even for the special cases where the background is purely 
magnetic or purely electric (ll], |l2| . 

However, the magnetic/electric cases are not the simplest ones which one can study in this 
context. As we will see in the following, the constant field background which leads to maximal 
simplification is the (Euclidean) self-dual one, given by 
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1 This Euclidean self-duality condition should not be confused with the similar condition a — ±fe which can 
be realized in Minkowski space and could also be called 'self- duality'. See [[L9| for a discussion of these various 
dualities. 
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This is because for such a field the square of the field strength tensor is proportional to the 
Lorentz identity (see eq. (|2.2| ) below) so that the inevitable breaking of the Lorentz invariance by 
the background field is in some sense minimized. In Minkowski space the self-duality condition 
requires either E or B to be complex. However, this does not imply that such backgrounds are 
devoid of physical meaning. Rather, they should be interpreted in terms of helicity projections 
|p0| . And indeed, it will be seen below that it is quite straightforward to extract, from the 
effective action for such a self-dual field, the low energy limit of the two-loop amplitude for the 
scattering of N photons with all equal helicities. In addition, we see at least four more good 
reasons for studying this particular self-dual case in detail: 

First, a detailed comparison at the one- loop level shows that the self-dual case with real B and 
complex E (corresponding to a real self-dual field strength /, where f 2 = \F^ V F^ V ', and called 
the 'magnetic' case in the following) leads to an effective Lagrangian whose properties are only 
marginally different from the ones of the magnetic Euler-Heisenberg Lagrangian, while the case 
with a real electric and complex magnetic field (corresponding to / imaginary, and called the 
'electric' case in the following) is a very good analogue of the electric case. This leads us to 
expect that the study of the self-dual case at higher loops may yield useful information on the 
generic properties of Euler-Heisenberg Lagrangians. 

Second, many quantities of physical interest in quantum electrodynamics are computable in 
Euclidean space. This includes the renormalization constants, and in particular the QED (3— 
function. 

Third, the constant self-dual Euclidean backgrounds that we consider here for QED generalize 
in a very simple way to the case of quasi-abelian self-dual constant backgrounds in QCD. Such 
backgrounds have been studied extensively in QCD as they have the special property that 
among the covariantly constant gluon backgrounds, only the self-dual quasi-abelian background 
is stable (at one-loop) under fluctuations. This has led to extensive studies of quark and gluon 



loops in such a background |21 



\, 23, 24]. Also at one-loop, much is known about QCD in 
the presence of arbitrary covariantly constant background fields. One-loop results along the 
lines of the Euler-Heisenberg formulas (1.1) and (1.2) can be found in [25, 26]. The purely 
chromomagnetic covariantly constant background is unstable [27, 21]. The renormalization 
of these one-loop results in pure QCD is complicated by infrared problems, problems which 
become even worse at the two-loop level [p^ [. 

A fourth motivation for studying constant quasi-abelian self-dual backgrounds is that they could 
provide useful information about effective Lagrangians in other self-dual backgrounds, in par- 



ticular instantons [3C, 31, 32]. 

As our main result in this paper, we will show that in the self-dual case, at two-loop, all parameter 
integrals can be done in closed form, yielding the following simple expressions for the two-loop 
scalar and spinor QED effective Lagrangians: 
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Here we have defined the convenient dimensionless parameter 
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where f 2 = jF^F^ ', as well as the important function 



£(x) 
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(1.7) 



with ^ being the digamma function ^(x) 
already presented in f33| . 



r'(x)/r(x). The scalar QED formula (|1.4|) was 



The extreme simplicity of these results, and the obvious similarity between the two-loop effective 
Lagrangians in the scalar and spinor cases, provides a new motivation for studying self-dual 
backgrounds. Our results provide a new example of the well-known connections among self- 



duality, helicity and super symmetry. Self-dual fields have definite helicity [20 J, and are closely 
related to supersymmetry (3(], |4|, 35 1. One consequence of this connection is that there exist 
remarkably simple formulas for loop amplitudes when the external fields have all (or almost all) 
helicities being equal. Such simplifications have been known for a long time in massless QCD at 
the tree level [p6f[ . More recently they have been studied extensively at the one- and two loop 
level, in abelian |37], 38, 39] as well as in nonabelian [4C, 41, 42] gauge theory. The close interplay 
been self-duality, supersymmetry and integrability has also been explored in this context. The 
simplicity of the structure of the QED/QCD amplitudes with all (almost all) helicities alike is 
thought to be deeply related to the integrability properties of self-dual Yang- Mills fields [^, f§4]] . 
In the present paper the connection between self-duality and helicity will be used to obtain, 
from the above results for the self-dual two-loop effective Lagrangians, the low energy limits of 
the QED N - photon amplitudes with all helicities alike. 

The outline of the paper is as follows. In section 2 we write down the one-loop self-dual Euler- 
Heisenberg Lagrangians for scalar and spinor QED, as well as for the 'quasi-abelian' special 
case in QCD. In section 3 we compute the two-loop corrections to these Lagrangians for scalar 



and spinor QED, using the 'string-inspired' worldline formalism [45, 46] along the lines of [J47], 
48, 15, 16 1 . This formalism is based on the representation of effective actions in terms of first- 
quantized particle path integrals, and has turned out to be highly convenient for computations 

H 0, 



49, 



51 



16 ) (see [53 1 for a review). In 



involving constant external fields [ 
section 4 we comment on the special properties of self-dual fields with respect to helicity and 
supersymmetry, and we use some of these properties to explain the similarity of our two-loop 
scalar and spinor QED results. We then proceed in section 5 to an explicit calculation of the 
low energy limits of the massive all '+' helicity N photon amplitudes, at one and two loops, in 
scalar and spinor QED. Section 6 contains a summary, and some possible future directions of 



work. In an accompanying paper [54] we study the weak- and strong-field expansions of these 
two-loop effective Lagrangians, and use them to test the techniques of Borel summation, and the 
associated analytic continuation properties and structure of the imaginary part of the two-loop 
self-dual Euler-Heisenberg Lagrangians for the 'electric' case. 



2 One loop self— dual Euler-Heisenberg Lagrangians 



2.1 Scalar QED 



For scalar QED, the (unrenormalized) one-loop effective Lagrangian in a constant F^ back- 
ground is 
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Here T is the proper-time variable for the loop scalar, and D is the spacetime dimension. We 
would like to evaluate this integral for the case of a self-dual field, F^ v = F^ w = \e ^y^F 01 ^ . 
The self-duality condition implies that 

F 2 = -fl (2.2) 
where f 2 = jF^F^, and 1L denotes the identity matrix in Lorentz space. Therefore 
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Renormalization involves subtracting the free field contribution, and also a charge renormaliza- 
tion, corresponding to subtracting the logarithmically divergent 0(/ 2 ) term. This leads to the 
following renormalized one-loop effective Lagrangian: 
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We remark that this integral can be expressed in terms of the function lnT(x) and its integral 
function: 
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where the function H(x) is defined as follows: 
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This one-loop result fl2.5|) may equivalently be expressed in terms of the Hurwitz zeta function 
(see, e.g., [55]]), but we prefer to use this representation here, since the two-loop result is a 
simple function of = S'(«). 



2.2 Spinor QED 



For spinor QED, the unrenormalized one- loop effective Lagrangian in a constant F^ v background 
is 



C^JF) = -2(4vr) 
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After renormalization, this leads to the following effective Lagrangian: 
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2.3 QCD 



In the nonabelian case, the natural notion of a constant field strength is covariant constancy, 
i.e. D^F^ = 0. This does not imply that F^ v is a constant matrix, so that further assumptions 
need to be made to arrive at Euler-Heisenberg type formulas. The simplest such case is the 
'quasi-abelian' one, where = A^T a , and therefore also F^ v , are assumed to point in a fixed 
direction in colour space: 



A%[x) = n a A,{x), F^(x)=n a F^(x) (2.10) 

The computation of the one-loop effective action with such a field is identical to the abelian one 
for the scalar and spinor loop cases, since the colour degree of freedom manifests itself only in 
a global colour trace. Thus for the self-dual case one obtains (compare (^) , (|2.7D ) 
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For the gluon loop one finds, in this quasi-abelian constant field case, the result |26|, 55, 15] 
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In these formulas, tr l denotes the Lorentz trace, tr c the colour trace, and we defined ip, = n a T a . 
For the self-dual case, ( 2.13] ) can be written as 
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(2.14) 



Let us specialize further to the case where all particles are in the adjoint representation and 
massless. We can then combine the contribution of m (real) scalars, n (Weyl) spinors, the gauge 
boson and its ghost (which gives minus the contribution of a complex scalar) into (putting 
D = 4) 
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For m = 6 and n = 4 we recover the well-known fact that the one-loop self-dual effective action 
vanishes for N = 4 SYM theory pi, (M M. 
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3 Two loop self— dual Euler-Heisenberg Lagrangians 



3.1 Scalar QED 



In [15, the two- loop Euler-Heisenberg Lagrangian in (Euclidean) Scalar QED was obtained 
in terms of the following fourfold parameter integral, 
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Here T and T represent the scalar and photon proper-times, and T a f, the endpoints of the photon 
insertion moving around the scalar loop. C a b and {y a ■ y&) are given by 
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They are expressed in terms of the worldline Green's function Gb an d its first an d second 
derivatives |l5| ]: 
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These are Lorentz matrices, and the above formulas should be understood as power series in the 
field strength tensor F^. The scalar function Gbi2 is given by Gbi2 = sign(ri— T2)— 2{t\— t<z)/T. 
Due to the translation invariance of those Green's functions one of the integrations / d,T a ^ is 
redundant, so that we will set t& = in the following. 

Removing the second derivative 5b by a partial integration with respect to r a or Tj, one can 
obtain the equivalent integral 
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We would again like to evaluate this integral for the case of a self-dual field. The worldline 
correlators (|3.3j) for such a field simplify to the following : 
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where Z^ u = eTF^ u . The first determinant factor is already known from the one- loop calculation, 
eq. (| 
scalar 



The Lorentz matrix C a b, being an even function of F, for this background becomes 
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Therefore in the self-dual case the T integrals in eqs,(|3.1|),(|3.4|) can be done trivially. In the 



same way as in the magnetic field case [15, 16], further simplification is achieved by taking the 
following linear combination of the representations (|3.1|) and (3.4), 
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(e = D — 4). We are only interested in terms of order 0(F A ) and higher in this Lagrangian, since 
the 0(F 2 ) term contributes only to photon wave function renormalization. For those terms the 
final T - integration is finite, so that in I\ we can set e = immediately. The same is not 
possible for I 2 , since this integral for D = 4 has a divergence at the points where the two photon 
end points become coincident, u = 0, 1. For the calculation of this integral we split it in the 
following way, 
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Note that the 8 a t, - term contained in G Bab has been deleted, since here the photon proper-time integral is a 
tadpole type integral which vanishes in dimensional regularization. 
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Of these four integrals the first one is of order 0(e) and can be omitted. The second and third 
ones are equal and elementary, 
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The fourth integral can be calculated by transforming to the variable 4> = (1 — 2u)Z, and using 
formula 22.10.11 of |58| (see [59] for a general technique for computing this type of integral). 
The result is 
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with £> the Euler Beta function. Expanding in e we find 
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We have written this result for the Lagrangian in terms of bare quantities, since it still requires 
renormalization. Mass renormalization requires us to subtract a term SmQ-^^C^f^ (fo), 

where £g2ff (fo) is the one-loop Euler-Heisenberg Lagrangian (|2.4|) , and drriQ the one-loop 
mass shift, both calculated in dimensional regularization. The mass shift is 
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we perform a partial integration on the first term on the right hand side of ( p. 13 ). After this 
the pole term in it takes the same form as the pole of the mass shift term above. Subtracting 
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the mass shift term we end up with the following expression for the renormalized Lagrangian 
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This can also be written as 
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Note that, since we have not yet performed the photon wave function renormalization, this 
formula is correct only beginning from 0(—^). 

For the calculation of this integral, we go back to the regularized version, eqs. ( |3.§| ) - ( |3.12| ), 
and note that it involves essentially only the integral 
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Including the mass renormalization term, and again disregarding 0(/ 2 ) terms, we can then write 
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Using ( |3 . 1 8| ) this yields, after an easy computation, our final result (|T 
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where we have finally subtracted out the terms of order 0(f ), 0(/ 2 ). 



3.2 Spinor QED 

The calculations for the spinor QED case are very similar, so we content ourselves primarily 
with a presentation of the results. 

For the two-loop Euler-Heisenberg Lagrangian in Euclidean spinor QED, the worldline formalism 



leads to the following parameter integral, which is completely analogous to fl3.4| ) flq , 16]: 
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The bosonic worldline Green's function Gb is given in (|3,3j), while the fermionic one is: 
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The calculation, together with the charge and mass renormalization, proceed in a manner anal- 
ogous to the two-loop scalar QED calculation in the previous section (for details of the mass 
renormalization, see |l5|l). After similar manipulations we find that the on-shell renormalized 
effective Lagrangian can be written as the following proper-time integral: 



C 



(2)(SD) 



m 



spin 



-a- 



167T 5 K Jo 



-2kZ 



dZ e 



+ ( log(2KZ) + 7 - - ) ( 1 



-Z + coth Z log 



sinh Z 



1 



1 d 
I^ZdZ 



_L d_ 

inZ dZ 



1 



coth z Z - — 7 - - + - + 



Z 2 



coth Z 
5 



6 24kZ 



(3.24) 



Here we have also subtracted the zero-field term and done the charge renormalization. 

As in the scalar case, this proper-time integral can be done in closed-form. After manipulations 



similar to those in the scalar case, we find the final result (|1.5| 



(2)(SD) 
spin 



(IC) 



m 



-a. 



(4tt) 3 k 2 



3e 2 ( K )-e'( K ) 



where £(/c) is the same function that was defined before in ( |1.7| ). This result (1.5) is remarkably 
similar to the two- loop scalar QED result (|0|) . We find it interesting that the spinor and scalar 
QED results for a constant selfdual background can each be written in such a simple form in 
terms of the same function This fact is discussed in the next section. 



4 Self-duality, helicity and supersymmetry at one- and two- 
loops 

In this section we explain why the renormalized scalar and spinor effective Lagrangians in a self- 
dual background are so similar, at both one- and two-loop. The self-duality of the background 
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has the effect that both the scalar and spinor cases are dramatically simplified from the results 
for a general constant background. But, on top of these simplifications, the scalar and spinor 
results ( |1.4| ) and ( |1.5[ ) are remarkably similar to one another, each being expressed as a simple 
function of The key to understanding this correspondence is the connection between self- 
duality, helicity and supersymmetry. It is well known that self-dual gauge fields are closely 



related to supersymmetry, and that they correspond to helicity eigenstates |3C| 3J, ^0, 35|. In 
this section we explain how these connections are manifest in our exact effective action results 
at both one- and two- loops. 

Consider the one-loop QED case. Since coth 2 t = 1 + l/sinh 2 £, we compare the one-loop 
renormalized effective Lagrangians (|2.4| ) and (|2.9| ) for the scalar and spinor cases, respectively, 
and find that 

/X">m--j££3 sd »<«) («> 

Note that this proportionality between the scalar and spinor effective Lagrangians is only true 
for a self-dual background. Furthermore, the relation (|4.1| ) holds for the renormalized effective 
Lagrangians, not for the unrenormalized ones. 

This proportionality between the one-loop scalar and spinor effective Lagrangians reflects a well- 



known supersymmetry of the self-dual background pO, 34 , EG, |3q| . The basic relation between 



self-dual fields and helicity can be traced to the following identity (we work in Euclidean space): 

W (i±2) = (D M (i + 2») (4.2) 

which holds provided the field strength is self-dual. This fact has the consequence pCj ] 
that the bosonic and fermionic operators 



A B = -D^D^ + m 2 (4.3) 

2 



A F = -(7 M D M ) 2 + m 2 = -D^D^- -o^F^ + m 2 (4.4) 



have identical spectra, but the fermionic case has a four-fold multiplicity (from the spinor degrees 
of freedom). This is due to a quantum mechanical supersymmetry p3] relating the bosonic and 



fermionic operators in (4.3) and (4.4). Therefore 



48L 8D) = \ log det A F = 2 log det A B = -2£^ SD) (4.5) 



which is precisely the one- loop relation ( |4.1[ ). 

Another equivalent way to see this is to note that because of the helicity/self-duality relation 
1.2|), the spinor propagator in a self-dual background can be expressed in a simple manner in 



terms of the scalar propagator in the same background [35]: 



S = - {ip- m) G (I±p) — GIp (j-p) + ^ (1 + JPGIP) (4-6) 
where the scalar and spinor propagators are 

G= -D 2 + m 2 ' S= lp+^l (4 ' 7) 
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This relation between the spinor and scalar propagators in a self-dual background provides 
another perspective on the relation (4.1). The spinor induced current can be related to the 
scalar induced current since, 

5 



spin/ 



/.(I) 

SA U spin 



= -ietri'fS) 

= 2ie{D li G + GD„) (4.8) 

where in the last step we have used (4.6). This current relation agrees with the one-loop effective 
Lagrangian relation ( |4.1| ). 

At the two-loop level, the supersymmetry relation ( |4.1| ) no longer holds, as can be seen by 
comparing the two-loop scalar and spinor results (|L4|) and (|1.5|) , respectively. Namely, 

e 2 (K) 



c 



-3m 



a 



(4.9) 



(4vr) 3 ,. 

This non-vanishing is because the two-loop effective Lagrangians are not simply logarithms of 
determinants. Nevertheless, it is interesting to see that each of the scalar and spinor two-loop 
effective Lagrangians can be written in such a simple form in terms of the same function 
To understand why this is the case, let us reconsider this computation in terms of more standard 
field theory techniques. Using Ritus's approach |L2| the two-loop effective Lagrangians can be 
written as [111 



C 



(2) 
spin 



(2) 
seal 



dx'V(x - x) tr [j^S(x, x')~f^S(x, x)] (4.10) 
dx'V(x - x') \d^G{x, x') D^G(x',x) + DpG(x, x') G(x', x)] (4.11) 



where S(x, x') and G(x, x') are the position space spinor and scalar propagators in the presence of 
the background field, and T>(x, x') denotes the internal photon propagator. We have not written 
tadpole terms, as these do not contribute after renormalization. These expressions dPcft and 
(4.11) are valid for a general background. Now, if we specialize to a self-dual background, then 
the spinor and scalar propagators become related as in ( |4.6| ), which means that the Dirac traces 
can be done in the spinor case, leading to an expression in terms of the scalar propagator: 

„2 



(2) 
spin 



dx'Vix 



-8DM(x, x') D^G(x', x) + l6D fl G(x, x') D„ G(x', x) (4.12) 



Thus, comparing with the scalar two-loop effective Lagrangian ( |4.1l[) , we see that the spinor 
two-loop effective Lagrangian for a self-dual background involves the same structures but with 
different coefficients. This is reflected directly in our explicit closed-form two-loop expressions 
(O) and (O). 



To make this connection more precise, we consider a general linear combination of these two 
structures: 



dx'V{x 



x 



A D^G(x, x') D^G(x', x) + B D^G(x, x') G(x', x) 



(4.13) 



where A and B are numerical coefficients. Now, in a constant self-dual background, the position 
space scalar propagator has the following simple form (up to a phase that is not important here) 



G(x, x ) 



ef 
(4vr) ; 



dt 



sinlr(i) 



cxp 



-2nt 



e f 



(4.14) 
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This means that we can express the general linear combination in ( 4.13| ) as 



m 



(4vt) 4 k 2 

4 

m 
(4tt) 4 k 2 



dt 



dt' 



ZdZe 



-2kZ 



e-Mt+t') j \{A + B) 4 coth (* + Q ~ j K 

sinh 2 (t) sinh 2 (t') \ (coth(t) + coth(i')) 2 ~ (coth(t) + coth(t')) 
'\{A + B) ^coth(Z) -Bk \ 



du 



o 



sinh 2 (Z) sinh(Z)(cosh(Z) - cosh(Z(l - 2u))) 



(4.15) 



where in the last step we have defined new integration variables Z = t + t', and u = ^hr. After 
renormalization, and using similar integration techniques as in section ||, the first term inside the 
curly parentheses in ( [4.15 ) leads to with a numerical coefficient \{A + B), while the second 
term leads to with a coefficient —(4 + Thus, the general linear combination in ( |4.13| ) 
leads, after renormalization, to 

J dx'V(x - x') [A D^G(x, x 1 ) D^G(x', x) + B D^G{x, x') B M G(x', x) 



m 1 
(4^)4^2 



{^ + B)e^)+ l -{A + B)i'(n) 



Now, the two-loop scalar QED effective Lagrangian ( 4.11 ) has numerical coefficients A 



(4.16) 

-B = l, 
On the 

other hand, the two-loop spinor QED effective Lagrangian ( 4.12 ) has A = —8 and B = 16, and 
an overall factor of which then leads immediately to our spinor result ( |1.5[) . Thus, the simple 
relation between the two-loop spinor and scalar effective Lagrangians for a self-dual background 



and an overall factor of — e , which then leads immediately to our scalar result ( |1.4| 



is ultimately due to the relation (4.6) between the spinor and scalar propagators in a self-dual 
background. 



5 The all '+' helicity amplitudes in the Euler-Heisenberg ap- 
proximation 

As is well known, the Euler-Heisenberg Lagrangian contains the information about the low- 
energy limit of the ./V - photon amplitudes at the given loop level. After specialization to a 
self-dual field it still contains information about the component of this amplitude which has 
all helicities alike. This component of the N - photon or gluon amplitudes generally leads to 



especially compact expressions Q36L 37, 38, p9, 4C, 41, [43, pOj. In this section, we will apply 



our results for the self-dual Euler-Heisenberg Lagrangians to an explicit computation of these 
amplitudes in the low-energy limit, at one and two loops. 

At the one loop level it is not difficult to perform a direct calculation of these amplitudes in the 
Bern-Kosower formalism. We present this calculation, too, since it is quite instructive to see 
how the equivalence between both approaches comes about. 



The Bern-Kosower master formula |45j, together with a judiciously chosen set of partial inte- 
grations |45|, |6^, pq| , allows one to write the one-loop N - photon amplitude in scalar QED in 
the following form [] : 

T^ l [h,e 1 ;...;k N ,e N ] = -e~ tm s s N ~ 2 / duj du 2 --- du N 



3 In this section we use Minkowski space (H ) conventions 



JO JO 
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xQ N (G Bi j) expi -is ^ G B ijki-kj\ (5.1) 

1 i<j=l ' 

Here Gsij =| Ui — Uj \ —(ui — Uj) 2 , Gsij = sign(uj — Uj) — 2{m —Uj). This representation holds 
on- and off-shell. It is important to note that the right-hand side represents the complete ampli- 
tude, not just a particular Feynman diagram. See [61, 53] for the relation of this representation 
to standard Feynman parameter integrals. The exponential factor exp{— is J2f<j=i Gsijki ■ kj} 
corresponds to the standard one- loop N - point denominator, while Qn corresponds to a nu- 
merator polynomial. This polynomial depends on the various Gg^'s as well as on the kinematic 
invariants Ei ■ kj, Ei ■ Ej. Qn can be decomposed in terms of 'cycles' and 'tails' as follows ^ 
Define the 'Lorentz cycle' Z n by 

Z 2 (ij) = -tr^) 
Z n {hi 2 ...i n ) = ti(F h F i2 ■ ■ ■ F tn ) (n>3) (5.2) 

where (Fi) uu = ki a Ei u — £iuki v is the field strength tensor for leg 'i'. Define the 'bi-cycle of length 
n\ G(hi 2 ■ ..i n ), by 

G(i\%2 ■ ■ ■ in) = GBi 1 i 2 GBi 2 i i ■ ■ ■ GBin^inGBin^Z nihil ■ ■ -in) (5.3) 
Then Qn has a decomposition of the form 

Qn = T,Qn Z2 "' Z ' (5-4) 

where the superscripts on a term indicate that it contains i bi-cycle factors of lengths z%, . . . ,Zi. 
Unless all indices 1, . . . , N in a given term are bound up in cycle factors there will be one more 
factor, called 'tail'. 

One advantage of this representation is that the polynomial Qn is homogeneous in the momenta; 
each term in it has exactly N factors of momentum. Therefore the low-energy limit of the am- 
plitude, corresponding to the case where m 2 becomes much larger than all kinematic invariants 
ki ■ kj, is simply obtained by replacing the exponential factor exp{— is J2i<j=i Gsijki ■ kj} by 
unity. The omission of this factor not only decouples the global proper-time integral, but also 
turns all terms involving tails into total derivatives. This leaves us with only those numerator 
terms which are products of bi-cycles. Thus the low energy, or Euler-Heisenberg ('EH'), limit 
of this amplitude can be written as 

r scai [ki,£i;...;k N ,£ N \ = — 16?r 2 m 2jV-4 — / dui du 2 '--J du N Q N '{G B ij) 

(5.5) 

where 

Q ( N H \GBij) = p ({*i*2 • • • i kl })P({ik!+i ■ ■ ■ ik 1+ k 2 }) • • • P({i kl +k 2 +...k j - 1 +i ■ ■ ■ iw}) 

part. 

(5.6) 

Here the sum is over all possible partitions of the set of indices 1,...,N into subsets, and 
P({iii2 ■ ■ ■ ik}) denotes the sum over all distinct bi-cycles which can be formed with the given 

4 We follow the notation of |63| rather than |53| . 
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subset of indices, e.g. P({ijkl}) = G(ijkl) + G(ijlk) + G(ikjl). The integral Jq 1 du\- ■ ■ dujq 
factorizes into copies of the basic 'chain integral' pjl 



/ duidu 2 ...du n GB12GB23---GBnl 

Jo 



in§n n even 



(5.7) 

n odd 



where B n denotes the nth Bernoulli number. The result of the integrations can therefore be 
written as 



/ dui [ du 2 -- - [ du N Q^ H \G B ij) 
Jo Jo Jo 



b^Fkxiih^ ■ ■ • 4i})frfc 2 -Ffc 2 ({4i+i • • • ^fci+fc 2 }) " " " bkjFkj ({*fci+fc 2 +...fcj-i+l • • • *Jv}) 

part. 

(5.8) 

where Fk({i\i2 ■ ■ ■ ik}) denotes the sum over all distinct Lorentz cycles which can be formed with 
a given subset of indices, e.g. Fi({ijkl}) = Z^ijkl) + Z^ijlk) + Z^ikjl). This allows us to 
write the amplitude in the following form: 



r (i)(EH) u , ] (-ie) N (N-l)\ /A v-^ j? rr- ■ ■ w 

F scal ! [kl,£i;...;k N ,e N \ = 167r 2 m 2jV-4 eX Pl 2^ ° 2m 1^ *2m({H*2 • • -12m}) 



m=l {ii...i 2m } 



(1...7V) 

(5.9) 



Here the second sum in the exponent runs over the set of all sets of 2m different positive numbers, 
and it is understood that after expansion of the exponential only those terms are kept which 
contain each index 1,2, ... ,N precisely once. 

It takes only a minor modification to generalize this computation to the spinor loop case: Ac- 
cording to the Bern-Kosower 'replacement' rules, the transition from scalar to spinor QED can 
be effected by replacing the above 'chain integral' (|5.7|) by the 'super chain integral' 



K — fn = J duidU2 ■ ■ ■ du n (Gb12Gb23 ■ ■ ■ GbvX — Gf12Gf23 ■ ■ ■ GFnl^j = (2 — 2") b n 

(5.10) 

As usual one must also multiply by a global factor of (—2) for statistics and degrees of freedom. 



Let us now compare with the calculation of the same amplitude via the effective action. In ( |1.2j ) 
we wrote the onedoop Euler-Heisenberg Lagrangian in terms of the two Maxwell invariants. 
Alternatively, the integrand can also be written in terms of traces of powers of the field strength 
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v m=l J 



According to standard field theory (see, e.g., p4J|) the low energy limit of the iV - photon ampli- 
tude can be obtained from the effective Lagrangian by replacing F with i J2iLi F%, expanding out, 
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and collecting all those terms which involve each F±, . . . , Fn, precisely once. The equivalence of 
( |5.9| ) and ( |5.11 ) therefore follows immediately from the combinatorial fact that 



tr 



{Ft + . . . + F N ) n 



all different r ■ ■ \ 



2n V F n {{hi 2 ...i n }) {2<n<N) 



(5.12) 



Here the sum runs over all possiblities of choosing a subset of n indices out of { 1 , . . . , N} . On 
the left hand side it is understood that terms containing the same F, twice are to be discarded. 

To obtain the on-shell helicity amplitudes, one now needs to construct polarisation vectors for 
all possible inequivalent helicity assignments, and compute the various e$ • kj's and e, • ej's. 
For small numbers of photons, this could be easily worked out for arbitrary assignments, using 



standard techniques such as the spinor helicity formalism |65|, gg, |67], 68, 36]. In this formalism, 



a polarisation vector with circular polarisation '±' for a photon with momentum k is written as 

4 = ± {q Z)^ k l ) (5.13) 



Here (q^ \ k T ) = u±(q)u^(k) etc. are basic spinor products, and q is a reference momentum 
which is different from k but arbitrary otherwise; see ||| for details and conventions. We will 
be concerned here only with the all '+' amplitudes (see, however, the last remark at the end of 
this section). For those, we know that all field strength tensors F\, . . . , Fn, and therefore also 
their sum F to t = Yli=i are self-dual tensors. Therefore, as in ( |2,2| ) the square of F to t must be 
proportional to the Lorentz identity. And indeed, using the spinor helicity formalism it is easy 
to see that 

Ff = 



where [ij] = (ki + \ kj ). It follows that 

-^tot = — /tot^ 

f tOt = Q ^ y [ 

l<i<j<N 



(5.14) 



(5.15) 



and 

tr[F t 2 ™] = 4(-ir(/ t 2 ot ) m (5.16) 

Thus, the all '+' component of the N - point amplitude can be obtained from the self-dual 
effective Lagrangian ( |2,4j ) as follows: Expand this Lagrangian in powers of /, i.e. in the weak 
field expansion; explicitly, 

Ai)(SD) m _ m A ^ (i) (2e/) 2 " 

^scal \J I — U^\2 c scal\'V ^An 



n=2 



J 1 ) y = B*n 

^ scal1 ' 2n(2n-2) 
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The N = 2n point amplitude is then obtained simply by taking the term involving / , replacing 
f 2 by / t 2 ot , and keeping only those terms in the expansion of f^ t involving each index precisely 



once: 



rW( Eff )[k £ +. -fcve+1 - ^ N c (1) fMv,u 

J-scal L^l) 5 ! > • • • ) rc AT,fcjvJ ~ //^2„2JV-4 C scalV"2VAiV 



= Jtot 
(f) 



(4-7r) 2 m 2 

all different 

{[12] 2 [34] 2 • • • [(TV — 1)7V] 2 + all permutations} 

(5.18) 



! 

2 ' 



N 

2T 



To generalize this result to the spinor QED and two-loop cases, all we need are the corresponding 



weak field expansions. For the one-loop amplitudes in spinor QED, we infer from (4.1) that they 
are proportional to the scalar QED ones: 

r«f H) [%,£+;... = -2r^ EH \k 1 ,et;...;k N ,e+] (5.19) 

At the two-loop level, we can easily obtain the weak field expansions of the Lagrangians ( |1.4| ). 
( |1.5| ) from the well-known asymptotic expansion of the digamma function at infinity (see part 
2). This yields 



rtWli, £ +. -kM£ + ] - air ^ c (2) f£Wy 



(4vr) 2 m 2 ^- 



(2) r x _ 1 / 2n - 3 | 3 ^4 B 2 k B 2n -2k 

seal W " (27r)2 \ 2 n-2 ° 2n ~ 2 2^ 2k (2n - 2fc) 



(5.20) 



for scalar QED, and 



r (2)( J E//) r , +. . fc +l _ g (ggT (2) m 

Jospin L ft l' £ 1 ) • • • J 'WV' e AN — Ztt7r ^4 7r p m 2Ar-4 C spinv~JAiV 



(2) , v 1 f 2n - 3 ^ S 2 fc #2n-2/c 

c s P inN " ^5 { 2n^2 + 3 2F (2n - 2fc) 



(5.21) 



for spinor QED. 



Let us discuss some implications of these results. In the massless case the one- loop all '+' N - 
photon amplitudes have been shown to vanish for N > 6 |57]]. Our result shows that this miracle 
does not generalize to the massive case at either one or two loops. 

The one-loop relation ( |5.19| ) is a direct consequence of the supersymmetry Ward identities for 
S-matrix elements first obtained by Grisaru, Pendleton and van Nieuwenhuizen |7(J. Those 
identities hold for massive as well as for massless supermultiplets in the loop. As was already 
clear from (4.9) this relation (|5.19| ) does not generalize to the two- loop level. The reason is, 



of course, that at two loops ordinary QED 'knows' already that it is not supersymmetric. The 
corresponding supersymmetry identity at the two-loop level involves, in addition to the scalar 
and fermion loop pieces ( ET20D , ( PT2TD , contribut ions from diagrams with a photino exchange 



as well as diagrams containing the scalar quartic coupling. At two loops our results could still 
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be used together with the supersymmetry Ward identity to infer the sum of those additional 
contributions to the 'all +' amplitude in the Euler-Heisenberg approximation. See the recent 
p9[ for an explicit verification of the vanishing of the two - loop four - photon 'all +' amplitude 
in massless N = 1 and N = 2 SUSY QED. 

Now that we have the above explicit all - N formulas, it would be of obvious interest to study 

(l 2) 

the large N behaviour of the coefficients spin (-f )• This will be done in detail in part 2. For 
now, let us just anticipate the facts that the leading asymptotic growth of the coefficients turns 
out to be factorial, and to be the same at one and two loops. 

Finally, let us mention that the same method of computation applied to the photon amplitudes 
with all but one helicity alike ( — !- + ■■■+) allows us to show without computation that those 
vanish in the low energy limit. Namely, gauge invariance implies that the Euler-Heisenberg La- 
grangian can, at arbitrary loop order and at an arbitrary fixed order in the weak field expansion, 
be written as a linear combination of terms involving only products of tr[.F n ]'s. After substi- 
tuting for F the sum of the individual F^s this yields sums of terms involving only products 
of Z n 's, where one of the factors must contain the Fi which has negative helicity, say, F\. But 
using the spinor helicity formalism it is easy to show that a tr (F^~ F^ ■ ■ ■ F^~) always vanishes 
(use ki as the reference momentum for all other legs). Therefore the Euler-Heisenberg limits of 
these amplitudes must vanish to all orders in perturbation theory. This argument applies, of 
course, to both scalar and spinor QED. 

6 Conclusions 

To summarize, we have computed the on-shell renormalized two-loop effective Lagrangian for 
both scalar and spinor QED, for the case of a constant self-dual Euclidean background. This 
case exhibits remarkable simplifications compared to the general constant background, and even 
compared to the constant magnetic (or electric) background. For the self-dual case the propa- 
gators simplify to such an extent that the final renormalized answers can be given in a simple 
closed form in terms of a common function which is essentially the digamma function. 

Although for (j) 4 theory an analogous result has been known for a long time jnj] we are not 
aware of any comparable two-loop result in gauge theory. 

Furthermore, the scalar and spinor QED cases are very similar to one another, a fact that 
we explained using the relation between self-duality, supersymmetry and helicity, which has 
the implication that the spinor propagator in a self-dual background has a simple expression 
in terms of the scalar propagator in the same background. At one-loop this implies that the 
scalar and spinor effective Lagrangians are proportional, but at two-loop the relation is more 
involved. This provides a new manifestation of the well-known simplicities arising in helicity (or 
self-duality) projections. 

Our explicit results for the self-dual two-loop effective Lagrangians have allowed us to compute 
the low energy limit of the massive N - photon amplitudes with all helicities alike, at one and 
two loops, for both scalar and spinor QED. This shows, in particular, that the vanishing theorem 
mentioned above does not extend to massive amplitudes. As a side result, we have shown that 
the amplitudes with all but one helicity alike vanish in the same approximation, to all orders in 
perturbation theory. 

From the mathematical point of view, we find it intriguing that, for both scalar and spinor QED 
and at both one and two loops, the final results for the effective Lagrangian turned out to be 
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expressible in terms of the 'mother' function H(re) and its derivatives. This property might be 
accidental, but we have reasons to suspect that, at least for QED in the quenched approximation, 
it will be found to extend to arbitrary loop orders. 

Moreover, we have seen that even in QCD the same would hold true at one loop if a gluon 
mass was introduced (say, as an infrared regulator). It would be very interesting to know the 
form of the result for the nonabelian case at two loops, i.e. the two loop effective Lagrangian 
for Yang-Mills theory in a constant self-dual quasi-abelian background. The two-loop worldline 
formulas derived in |29] for the general nonabelian Euler-Heisenberg type Lagrangian should 
provide a good starting point for its computation. 

(2) 

Our results also suggest a new approach to the computation of the Lagrangians <£ sca j S pj n fo^" 
a general constant background. Rather than performing a double expansion in the variables 
a, b, as has been done in [12, pj]], it might be preferable to use an expansion in the variables 



H = a + ib and H = a — ib. The integration techniques which we have used here for the self-dual 
case, corresponding to H = 0, would in principle also allow one to obtain, at any fixed order of 
this expansion in H n , a closed form expression for the coefficient as a function of H. This might 
lead to progress even in the computation of the magnetic and electric cases. 

We anticipate that the techniques developed here may also help to extend the one-loop instanton 
background calculations (3^, 31, ^] to the two- loop level. 



In part 2 of this paper [^J], we investigate the structure of the imaginary part of the two- 
loop self-dual Euler-Heisenberg Lagrangians for the 'electric' case. We also study the weak- 
and strong- field expansions of these two- loop effective Lagrangians, and use them to test the 
techniques of Borel summation, and the associated analytic continuation properties. 
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